An Expansion for Laplacian Integrals in Terms of 

Incomplete Gamma Functions, and 

Some Applications* 

By EDWARD C. MOLINA 

Introduction 

LAPLACE has given us, in the Theorie Analytique des Probabilities, 
Book I, Part II, Chapter I, a method of approximating by means 
of series to the value of a definite integral of the type 



1= PyiV* •••?»'•***. 



where yt, yi, • • ■ y n are functions of x whose exponents 0\, 02, • • • 6 n 
are of the same order of magnitude as a large number 6. The last 
function in the integrand, #, embraces all factors whose exponents 
are of low order of magnitude compared with 6. The integral here 
considered must not be confused with the well known "Laplacian 
Transform" integral which is also embodied in the Theorie Analytique. 
When the function $(x) is other than a mere constant, the Laplacian 
method as presented in the Theorie Analytique does not, in certain 
cases, give us a series which reduces to its first term as approaches 
infinity. The object of this paper is to present a modification of the 
method which gives the desired result and to present two applications 
of considerable practical importance. The modification consists in 
divorcing the function from the factors of the integrand raised to 
high powers and associating <f> with the factor (dx/dt) which makes its 
appearance with the Laplacian change of variable from x to /. 

Deduction of Modified Expansion 
Setting 6jd = r, we have 

yfiyft . • • y n g n = (yfi^'i . . . y B '«)*= /, 

(1) 1= r^(x)J<t>(x)dx, 

which we shall write in the form 



1 - P !>(*)"*.(*)<**.] 



* Presented by title at International Congress of Mathematicians, Zurich, Swit- 
zerland, September, 1932. 
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where N = 6 + w, <£„,(*) = [y(x)]~"^(«). Introducing the parameter 
w does not constitute an essential modification of the procedure 
given by Laplace, but we shall find that w plays an important part 
when we come to applications of the expansion presented in this 
paper. 

In what follows it is assumed that y(x) is a positive monotonically 
increasing or decreasing function in the range x h x 2 ; the extension to 
a range of integration divisible into subranges within which y(x) is 
monotonic will be obvious. 

Let X be that one of the two limits x u x 2 for which y(x) has its 
greatest value; set y(X) = Y. Assume with Laplace that 



1 



(2) [log Y - log yix)^ 1 = (x - X)g(x, X), 

u being a positive number or zero and g a function of which (x — X) 
is not a factor. Set 

(3) y = Ye-' u+1 , 

t t 2 t 3 

(4) <t> w (x) (dx/dt) = A + A x p + A 2 ^ + A 3 ^ + 

These two equations give (certain well-known conditions being 
fulfilled) 

(s) i=f ^-j;-^- 

Finally set Nt u+1 = T and we obtain the modified Laplacian expansion 

(6) / = r*T o (V) H 'BstP(n„ r 2 ) - P{Ih, Tijf} 

wherein 

Ha= i±l t 

r, = N[)og Y - log y(«0], =0 K * = *tf 
T 2 = iVQog F - log y(Xi)2, =0 if X = x 2 ; 

T H ~ l e- T dT 

V = Tx or T 2 . 
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For the computation of P(H a , T') recourse may be had to the 
extensive "Incomplete Gamma Function" Tables edited by Karl 
Pearson which give the ratio of the incomplete to the complete function. 
When //, is a positive integer we have 

* 7-/*_-r/ 

P(ii„T') = Y. L 4 r - t 

the well known Poisson Exponential Binomial Limit for which short 
tables will be found in Pearson's "Tables For Statisticians And 
Biometricians " and in T. C. Fry's "Probability and Its Engineering 
Uses." 

Applications 
I 
Consider the incomplete Beta Function 



Jo 



for positive integral values of n and c such that n is large compared 
with c. Its modified Laplacian expansion is 

I(P) - (P - 1)!(» + w)- c Z(n + w)-">A(m, c - l)P(c + m, T t ), 

ro=0 

where T 2 = — (n + w) log (1 — p) and, setting w = (d — c -\- 1) /2, 1 
A(0,c- 1)= 1, 
4(l,c-l)«-(^, 

.(2, c - 1 ) = ( c t 1 )[^4±-], 

^ g -i)-( g + 2 )[ (c -y + 4 

i4(4, c- 1) = ( C + 3 ) \ 5(C - 1)2 ~ (2 ~ 3 ° d2)(c ~ 1} + l5di ] 

As many more coefficients as one desires can be obtained by means 
of the equations 

A(m, c-l)= ±Q(m, k) ( ° " * + "'). 
t =o \ w + « / 

(?(ro + 1, k) = (w + k)Q(m, k - 1) + (w + Jfe)Q(m, fe), 

G(*. o) - 1. 

1 See Appendix II for the details of the analyses. 
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The expansion given above for I(p) is of great practical value in 
connection with the evaluation of the incomplete binomial sum- 
mation. We know 2 that 

P(c;n,p) = £Qp*(l -PY~ X 
therefore, we have 



n\ ^ / 1 



This approximation formula is submitted for ranges of values of c 
and n such that c is of the order of magnitude of yjti or less. But for 
values of P{c\ n, p) which are of most interest, say, 

.0001 < P < .9999, 

c is of the same order of magnitude as 

a — np 

when n is a large number. Therefore, the formula is particularly 
valuable for values of p which are of the same order of magnitude as 
1 f-yjn or less. 

The first table at the end of this paper indicates the degree of 
accuracy obtainable by taking d = and using only one or two 
terms of the approximation for P(c; n, p). The table also indicates 
the result of taking d such that T 2 = a = np, and using only one, 
two, three or four terms of the approximating series. 

1 am indebted to Miss Elizabeth McCusker and Miss Catherine 
Lennon of the Department of Development and Research of the 
American Telephone and Telegraph Company for the computation 
work involved in the construction of the tables at the end of this 
paper. 

II 

Consider the Bessel Function 

J»0) = ,-lf } 1 n r ( sinx ) 2n eMC08X) dx > 

V «T(» + 2) Jo 

2 See Laplace, " Theorie Analytique des Probabilites," 1st edition, p. 151. 
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for b large compared with n, n being a positive integer. Writing 



e-"I n (b) = -= 



V?r£ + 1) f <» *V*r+*~r*« 



leads us to the modified Laplacian expansion 

2b 



\7r(26 + w) 
where 

1 \" i4i»+8» P(» -I- m + |, 26 + w) n\ 



=o\86 + 4wJ 2 (2») ! P(» + i 26 + w) (» + m) ! 

To determine the coefficients -4 2n +2m we proceed as follows: 
For the integral now under consideration 

N = 2b + w, *i = 0, X2 = it, 

y( X ) = g-(i-*«W) j X = 0, F= 1, 

(log F - log y)i - (| — | cosx)i = *f (*), 

gf» = - sin - = £ 



x 2 t f 2 2 * +l (2£ + 1)!' 
so that 

u= I, y(x) = e~ l \ T x = 0, T 2 = 2b + w, 

0»(*) = (i sinx) 2 "e w <*»>', 

B 2n+2m = [n!r(» + |)/(w + m)!(2fi)!2*^i]i4 1 » K ,. 

The relations above give 

<t> w (x)(dx/dt) = 2/ 2 »(l - t 2 ) n -h»'\ 
Let 

XM' 2n W 2 " +2 "/(2" + 25)! = 2/ 2 "(l - **)»-*, 



so that 



A' 2n = 2(2»)!, 

2rc - 1)(2» - 3) 

2'-'s! 



,, _ (- 1)'(2». - 1)(2» - 3) • • • (2» - 25 + 1)(2» + 25)! 

■^ 2»+2« — 



then 



4> w {x)(dxldt) = e^-ZA'w, 



,t&" *"-*•" (2» + 25)! 

E^ n+2i+2 ^'2n+2 S wV^! 
s = 

3 An expansion of this type for Bessel functions was derived many years ago by 
Hadamard; see Watson, "A Treatise on the Theory of Bessel Functions," pp. 204 



E E/ 2 " +2s+2 M' 2n+2s ^-/*!(2w + 25)! 

;.-=o s=o 



and 205. 



— ).P(3/2, 26 + .5)/Vtt(26 + .5)1 Si 
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and 

m 

A 2n+2m /(2n + 2m)\ = £(w k IW)A' 2n+2m - 2k l(2n + 2m - 2k)\ 

We must now assign a value to w. One which suggests itself is to 
take w such that the second term in our expansion for 7„(6) vanishes; 
in other words, such that A 2n+2 shall be zero. This merely requires 
that we set w = n — 1/2, giving 

A 2 „ = 2(2»)! f A 2n+2 = 

A 2n+ i = - (2n - 1)(2« + 4)!/2, A 2n+6 = - (2*1 - l)(2n + 6)!/3 

A 2n+8 = (2n - l)(2?i - 5)(2» + 8)!/16 f etc. 

For n = 1 and w = 1 /2 we have 

A 2 = 4, ^4=0, A 6 = - 360, ^ 8 = -- 13440, • • • 

r*m = [(2 2 * 

wherein 

5, - 1 - 15(86 + 2) p(3/2( 26 + .5) 

- 140(86 + 2) p(3/2) 26 + 5) 

The degree of accuracy obtainable when only one, two or three 
terms of Si are made use of is indicated by the second table at the 
end of this paper. For some values of 6, ranging from 6 to 16, com- 
parison is made between the successive approximations for e~ b Ii{b) 
and its true value as given by Watson in his Theory of Bessel Functions. 
The last three columns of Table II give the figures for w = instead 
of w= 1/2. 

Appendix I 
When the expansion of 4> w (x)(dx/dt) is not obvious we may proceed 
as does Laplace in expanding (dx/dt). 
Equations (2) and (3) give 

(7) x = X + tg~\ 

This last equation, together with the Lagrange-Laplace expansion 
theorem, gives for a function of x the expansion 



s=l J - 



d 
dx 



S(i)'™L- 
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But df(x)/dt = f(x)(dx/dt). Therefore, taking /(*) such that /'(*) 
= <t>w(x), 

. , . dx -s-^t' 



or 



(<s) 



A . = [D/r^'iWk, D *' = dk" 

= y- w zK[d x >(x - X)Wr (l+1,+W) '^)U 

t=0 "! 

4f« y- V w ' [/)x'- ( " +1 > p 0(^)g- ( ' +1)+( " +1)p l = , Y 



S! 



.-o^ 



[5 - (« + 1>]! 



.Se/a Function 
For the incomplete Beta Function we have, 

X=0, u=0, Y=l, g(0)=l, 



Therefore, 
A. 



a— (e— 1) w « 



*l ,?0 *l£ 



c - 1 + r\ [D x '~ v x 



»— » Y c— l+r a — (s+l— 



•>l-o 



(*-*)l 



or 



Aj = -WwV-lC-p-V t C _ 1 _]_ f \ [- £ ) il l-( C -l)-V-rg-(,+l-,)-| x=0 

5! = £i v! r to \ r ).\s - (c - 1) - v - rll * 

Since A a = for s < c — 1, set 5 = c — 1 + m/, then interchanging 
z; and m — v we obtain 



(9) 



A c . 



c—l+m 



m niim—V 



= E 



vr 



c- 1 +r\ [P x — ■g- (e+,) l-o 



(c - 1 + w») ! ,f (/» -f)!r=o\ c — 1 / (w — r) ! 
To evaluate (9) we may have recourse to the formula 4 



taking 



(-l)*( J ]£)zW 
(5 + K)g*+* 



M = v — r, 

S = c + v. 

4 See formula 25, page 15, of Schlomilchs Hoheren Analysis, Band II, Dritte 
Auflage, 1879. I am indebted to Mr. J. Riordan of the Department of Development 
and Research, American Telephone and Telegraph Company, for calling my attention 
to this formula. 
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Bess el Function 
For the Bessel integral with b large compared with n we have 

X = 0, 21= I, g(0) = 1/2, 7=1, 

0(.r) = [[sin (§)]*[*- sin2 (f )]}"= K**) 1 - ta) 4 3" 

= Z o (^)(-i) r W 2n+2r - 

Substituting in (8) we have 

il. _W* nr /w\ [D J '- 2 ".r 2 »+ 2 -g- (, ~ i! "~ 2 ''~ 2r+1) ]x = o 

where 2.R > s - In - 2v\ note that if R > n, then ( n ) = 0. 

Since A „ = for s < 2w and g is an even function of x, set s = 2n 
+ 2w. Then Leibnitz's theorem and interchanging v with m — v give 

(11; (2 W + 2m)! .toCm-*)!.^ ' W (2»-2r)l 

where again we may have recourse to (10) for the differentiations of 
negative powers of g(x). 

Appendix II 
Writing the incomplete Beta Function in the form 

I(p) = f P (l - .%-) n+u, * c-1 (l - .r)- (c+u,) </x, 
Jo 

we now have 

N = n + w, xi = 0, xi = p, 
y(x) = (l-x), 7=1, X=0, 
<f> w (x) = x'-^l - x)~ ic+w \ 
log Y -logy = - log (1 - x), 

= •%-£(*), *W = ( i +5 + J "•). 
u = 0, 
y{x) = e~ l , 
T,= 0, 
T 2 = - (n + w) log (1 - p). 
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From these equations we derive 
<l> w (x){dxldt) = e"'(e' - \) c ~ l 

-£(a)[I(-» , ( e I 1 )<-+«- 1 -H 1 

s — c — 1 + m. 

Therefore (see any work on finite differences), 

A e - X + m = b«-l( W )o-l+m = (_ l)m A c-i(_ w _ c _|_ 1)0-1+*. 

Now 

A c-l( w )c-l+ m +l - , t ,A« ; -l(w) c -l+'» + (c - 1)A C - 2 (W + l)c-l+m 

= (w + c - l)A c - I (w) e - 1+m + (c - l)A c - 2 (w) c - 2+m+1 . 

Increasing c by 1, decreasing w by 1 and setting 

A c _ 1+m = (c - l)\A(m, c - 1), 
we have 

A(m, c) — A{m, c — 1) = (w + c)^4(w — 1, c) 
or 

Aei4(w, c — 1) = (w + c).4(m — 1, c) 

where the subscript c implies that A c operates on c. Assume that 
A(m, c - I) = ± Q(m, k) ( C ~ l . + , m ) , 

where the coefficients Q(m, k) are functions of w; then 

c-\-m\ 
fe^b w ' ' ' \ W 

=]|W)[^ 

='2 I [(w+^)<2(m, *-i)+(w+A)e(«, *)]( c+ ™ 



A c /l(m + l, c-l)= ZQ(m, k)(w+c) \ ^ +jfe i 



t=o _ \w+^ 

from which, by finite integration with reference to c, we obtain 

^(«+l,c-l)=^ 1 [(w+ife)g(w,ife-l)+(a;+ife)(?(w f *)V C ~V , ;^t 1> ) 

t=o \ ?n-\-l-\-k I 

so that 

Q{m+ 1,0) = W (?(m, 0), 

<?(w + 1, *) - (w + *)(?(/», k - 1) + (w + k)Q(m, k), 

< k < m + 1, 

(?(m + 1, w + 1) = (2m + l)(?(w, m). 

Since 

/i(o, c -i) = 1 = ( c o 1 ). 
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we know that 0(0, 0) = 1. Hence 

4(l,c-l) = «;( 1 C ) + ( 2 C ). 

^a-D-^t'J + ^ + ^CtO + ^t 1 ) 



c + 2 



4(3, c - 1) = w 8 ( " 3 " ) + (6w 2 + 4«> + 1) v , 



c + 2 



+ 15, ! -to ( C t 2 ) + 15 ( C 6 2 



4(4, c - 1) = w*( C ^ 3 ) + (10w» + 10w 2 + 5w + 1) (' | 3 ) 

+ (45w 2 + 60w + 25)( C ^ 3 ) 

+ (lOSw + 105) ( C + 3 ) + 105 ( C + 3 )- 
etc. 

As yet the value of w has not been specified. Since when m is an 
odd number and w = - (c - l)/2, the equation 



gives us 



Ac-l( w )c-l+m _ (_ 1)M^«-1(— W — C + 1) 

/ r — 1 \ c - 1 + m 
A-' I - ^^ I = 0, 



e— 1+m 



let us write w = (d/2) - (c - l)/2, leaving d arbitrary; we obtain 



4(0, c - 
4(1, c- 

4(2, c- 

4(3, c - 

4 (4, c - 



= 1, 



1J2' 

c+ 1\ [(p ■ i: -I- .vr 
2 



)[- o 

c + 2 \ f (c - \)d 4- <f 3 



] 



)[ 



c + 3\ r 5(c - D 2 - (2 - 30<j 2 )(c - 1) + 15(^ 1 



4 / I - -i 11 

etc. Finally, for d = 0, or w = — (c — l)/2, 

4(0, c- 1) = 1, 4(1, c - 1) = 0, 

A(2,c-l) = ( C + l )l> 4(3, C -l) = 0, 

4(4,<-l) = ( J )(-jg- 
etc. 



4(5, c - 1) = 
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TABLE II 

Bessel Function 
e~ b h{b) 



v> =* — 1/2 - 1/2 for n = 1 



(***) 



'P -Pi^ 



(^) 



PP) 



<P -Pi\ 



10 
12 
14 
16 



.152052 
.134143 
.121263 
.111464 
.103698 
.097350 



.007497 
.004043 
.002529 
.001731 
.001259 
.000958 



.001457 
.000587 
.000293 
.000167 
.000104 
.000070 



-.000331 
-.000098 
-.000038 
-.000018 
-.000010 
-.000006 



.071134 
.051474 
.040358 
.033198 
.028198 
.024510 



-.004188 
-.002187 
-.001344 
-.000911 
-.000658 
- .000497 



-.000701 
-.000261 
-.000125 
-.000070 
- .000042 
-.000029 



P = value of e b Ii{b) from Table II of Watson's Theory of Bessel Functions. 
/ p _ p \ 
I p — - J = proportional error incurred by using the first m terms of the series 

for e~ b I\{b) presented in this paper. 



